Size dependence of the internal energy in Ising and vector spin glasses 
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We study numerically the scaling correction to the internal energy per spin as a function of system 
size and temperature in a variety of Ising and vector spin glasses. From a standard scaling analysis we 
estimate the effective size correction exponent x at each temperature. For each system with a finite 
ordering temperature, as temperature is increased from zero, x initially decreases regularly until 
it goes through a minimum at a temperature close to the critical temperature, and then increases 
strongly. The behavior of the exponent x at and below the critical temperature is more complex than 
suggested by the model for the size correction that relates x to the domain- wall stiffness exponent. 
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Since Gibbs, thermodynamic transitions have been 
classified according to the critical behavior of the spe- 
cific heat, or equivalently of the critical temperature de- 
pendence of the internal energy. One of the disturbing 
features of the spin-glass transition has always been that 
there appears to be no thermodynamic signature whatso- 
ever of the critical temperature, except in the mean-field 
limit. In finite dimensions, the specific- heat exponent 
a is strongly negative and the energy changes perfectly 
smoothly as a function of temperature through the crit- 
ical temperature T C X 

The standard scaling form for the finite-size correction 
to the internal energy per spin e(L) is 



e(L) 



aL 



(1) 



where L represents the system size. An exponent #e can 
be defined by 9e = d — x (d represents the space dimen- 
sion). In Ising spin glasses (ISG) it has been surmised 
from "droplet" arguments 2 *^ that this correction is di- 
rectly related to the energy associated with domain walls, 
for which independent numerical measurements can also 
be carried out. Thus it is expected that at zero temper- 
ature 9e is identical to #d\¥; the domain- wall stiffness 
exponent. This conjecture is related to the controver- 
sial question of the form of the elementary excitations in 
spin glasses, and the identity should be valid if periodic 
boundary conditions simply introduce supplementary do- 
main walls. It is exact for Migdal-Kadanoff spin glasses*^ 
While the argument was introduced for zero tempera- 
ture, it has also been invoked for finite T£ At a con- 
tinuous transition the singular part of the free energy 
divided by the temperature scales as length~ rf . Because 
(T - T c ) ~ length" 17 " (see Refs. H and@), if T c > 0, 



x(T = T c ) = d- 1/v 



(2) 



If T c = 0, 6> DW = -l/v at T = and x(0) = d+l/v£ 

We have carried out Monte Carlo measurements of the 
size dependence of the energy as a function of tempera- 
ture in the mean-field Sherrington-Kirkpatrick (SK) ISG 
model, in the Edwards-Anderson ISG with Gaussian in- 
teractions in dimensions 2, 3, and 4, in the gauge glass 



(GG) in dimensions 2, 3, and 4, and in the XY spin 
glass (XySG) with Gaussian interactions in dimension 
4. In all systems with a nonzero ordering temperature, 
x(T) is strongly temperature dependent below as well 
as above T c . The effective exponent initially decreases 
progressively as T increases from zero; it passes through 
a minimum at a temperature T m j n close to T c , and from 
then on increases sharply. The data for the mean-field SK 
spin glass and for the finite-dimensional systems follow 
strikingly similar patterns. We associate the observed 
minimum with the critical behavior in Eq. (j2J which po- 
tentially provides a powerful criterion for identifying T c 
in spin glasses from purely energetic measurements. 

The value of the stiffness exponent in the d = 3 gauge 
glass has been source of controversy: results have clus- 
tered either close to 6 >nw ~ 0.0 (Re fs. HHTol) or close 
to 6» DW « 0.27 (Refs. 



II 1 111 12111 311 14T ). This situation has 
been analyzed by Akino and Kosterlitzjii*!^ who show 
that following the boundary conditions imposed, either a 
"best twist" (BT) value (near 0.27) or a "random twist" 
(RT) value (near zero) of #dw is obtained. They asso- 
ciate the BT value with domain walls, but they say "We 
do not understand what, if anything, 9b,t means...," im- 
plying that this 9 could be nonphysical and simply an 
artifact arising from an inappropriate choice of bound- 
ary conditions. If interpreted in terms of a domain- wall 
stiffness exponent, our T — estimate is compatible with 
6>dw ~ 0. 

In the canonical [Ising, XY] Edwards- Anderson spin 
glass [Ising, XY] spins on a hypercubic lattice of size L 
interact with their nearest neighbors through random in- 
teractions whose strengths follow a Gaussian distribution 
with zero mean and standard deviation unitji 



H 



Jij^iSj 



(3) 



Periodic boundary conditions are applied. The mean- 
field limit system is the SK model. In the gauge glass, 15 
XY spins on a hypercubic lattice of size L interact 



2 



through the Hamiltonian 



H = —j'^j cos( 



- Aij), 



(4) 



where the sum ranges over nearest neighbors. <pi repre- 
sent the angles of the spins and Ay are quenched ran- 
dom variables (gauge fields) uniformly distributed be- 
tween [0, 2ir] with the constraint that Ay — —Aji. J = 1 
and periodic boundary conditions are applied. 

In the present series of simulations, samples 
are equilibrated using the parallel tempering Monte 
Carlo method-iSJi The sizes studied are N — 
36,64,100,121,144, and 196 for the SK model, L = 3 

- 6 in the d = 4 ISG, L = 2 - 5 in the d = 4 GG, L = 3 

- 5 in the d = 4 XYSG, L = 3 - 6, 8 in the d = 3 ISG, 
L = 2 - 6, 8 in the d = 3 GG, L = 3, 4, 6, 8, 10, and 12 
in the d = 2 ISG, and L = 6, 8, 12, and 16 in the d = 2 
GG. 

The equilibrium energies e(L,T) for a given system 
are averaged over at least 10 3 disorder realizations for 
the largest system sizes. Due to small differences in the 
random interactions, there are sample to sample fluctu- 
ations at each size. At each temperature, the data for 
each system are fitted using Eq. Q), with e 00 (T'), a(T), 
and x(T) as temperature-dependent fitting parameters. 
Note that for each value of L, exactly the same set of 
samples is studied over the whole range of temperatures. 

We first discuss the mean-field SK model and then the 
finite-dimension models in order of decreasing dimension. 
In the SK model, x(T) — 2/3 at the ordering temperature 
T c = l, 18 in agreement with Eq. ©. For T > T c the 
power law [see Eq. lfl]l] is an approximation to a sum 
of terms in N~ k with k E N. We fit e(N,T) assuming 
that the power law with exponent x(T) is valid at all 
temperatures below T c = 1 and is a useful approximation 
above. The infinite-size limit energies eoo(T) have been 
established numerically to high precision^ and we use 
these values as input in our fits. In Fig.^ it can be seen 
that x(T) is temperature dependent and passes through 
a minimum at T « 0.85. At T = 0, x(T) again tends 
to a value very close to 2/3, consistent with a directly 
measured zero-temperature estimated 

For the d = 4 ISG (Fig. HJ, E (O) = 0.71 ± 0.08 , 
in agreement with estimates for the T — domain-wall 
stiffness of the bimodal d = 4 ISG: 6» D w = 0.65 ± 0.04 
(Ref. i3) and 0.82 ± 0.06 (Ref. However, as the 

two estimates are fairly different from each other, a more 
stringent comparison must await a definitive estimate for 
#dw(0)- There is a minimum in x(T) at T ~ 1.5, a 
temperature lower than T c = 1.80 ± 0.02 (Refs. |2^ and 
H). We find x(T c ) = 3.10 ± 0.05, which is higher than 
d- \jv = 2.8 ± 0.1. 

For the GG in dimension 4, our data only extend down 
to T = 0.7; x(T) has a minimum at T « 0.85, slightly be- 
low the ordering temperature (T c = 0.89 ±0.01, Ref-HJ). 
By extrapolation, we estimate 6*e(0) = 0.54 ± 0.05 and 
x(T c ) = 3.42 ± 0.02, which, with v = 0.70 ± 0.1^ is 
distinctly larger than d — \jv= 2.6 ± 0.2. 



0.9 - 



0.7 



0.6 



0.5 



SK 

T min « 0.85 



■ 



0.4 



0.8 



1.2 



FIG. 1: Data for x(T) for the SK model. The data show a 
minimum at Tmin ~ 0.85 and x extrapolates to the mean-field 
value for T -» 0, 2/3. The point at T = T c = 1 is the exact 
result (Ref. El). In this and all subsequent figures the dashed 
line is a guide to the eye and the arrow marks the position of 
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FIG. 2: Data for x(T) for the d = 4 ISG. The data show 
a minimum at T m i n ~ 1.5, close to standard estimates of T c . 
We find E (T -f 0) = 0.71 ± 0.08 . 



We have also analyzed data on the four-dimensional 
XYSG with Gaussian interactions (see Fig- - The cor- 
rection exponent x(T) behaves similar to the other d = 4 
cases, with 6>e(0) = 0.60 ± 0.05 and a well defined min- 
imum at T m j n = 0.67 ± 0.02. We are not aware of mea- 
surements of #dw(0) or of T c for this system, although 
T c « 0.95 has been reported for the four-dimensional 
AY"SG with bimodal interactions^ 

Figure 0] shows a similar behavior for the d = 3 
ISG. Again, x(T) initially decreases as temperature in- 
creases. The data are consistent with two tempera- 
ture points reported for T = 0.7 and T — 0.8 by Ko- 
mori et al£. The extrapolated low-temperature limiting 
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FIG. 3: Data for x(T) for the four- dimensional XYSG. As we 
have the same number of system sizes as fitting parameters, 
we cannot estimate an error bar to the fits. The data show a 
clear minimum at T m ; n = 0.67 ± 0.02. 



1 1 1 [ 1 1 


1 ■ 


XY (d = 4) T 


w 0.67 / ' 

min . 


- 
■ 


■ 


■ 




■ 


■ 


1 " \ 


/ / j 




■ 


•• 

i 1 i i i i i 


iii 


0.25 0.5 


0.75 1 


T 





3.6 



3.4 



3.2 - 



GG (d = 3) T . » 0.45 

v ' min -,- 1 




i 



0.3 0.4 0.5 0.6J 

».„..•■'" 

I I I I I I I I I I I I I I I L 

0.2 0.4 0.6 0.8 



FIG. 5: Data for the three-dimensional gauge glass x(T). The 
data show a minimum at T mln « 0.45. In addition, x(0) fa 
2.99. In the inset the data are enlarged around the minimum. 



3.4 




FIG. 4: Data for x(T) for the Ising spin glass in dimension 3. 
We see a minimum which coincides with T c [this is enlarged 
in detail in the inset of the figure], T m i n » 0.89, kb(0) « 2.85. 



value is 0e(O) = 0.15 ± 0.02. This is consistent with a 
zero-temperature measurement 6e{0) = 0.135 ± 0.037. 4 
The values from the three independent studies taken to- 
gether suggest a zero-temperature limiting value E (O) 
that is significantly lower than the directly measured 
three-dimensional domain- wall exponent value 0dw(O) = 
0.19 ± O.Oli 20 ! 27 There is a minimum in x(T) at a tem- 
perature T m j n « 0.89 which, in this case, agrees within 
the quoted errors to T c estimated independently by other 
techniques,^ 2 - 2 ^ T c = 0.94 ± 0.03. Using v = 1.65 ± 0.1 
(Refs. I|3 and|H, d-l/v = 2.4 ± 0.05, a value below 
x(T c ) « 2.75. 

Figure [S] shows x vs T for the d — 3 GG. Here, the 
low-temperature limit corresponds to 0e = 0.010(12). 



As in the other cases, x(T) decreases with increasing T. 
There is a minimum in x(T) at T w 0.45, which again, 
within the error bars, agrees with the ordering tem- 
perature T c estimated from other methods We find 
d-l/v = 2.28 ±0.03 {v = 1.39 ±0.03, Ref. HI), whereas 
the measured value is much larger, x(T c ) = 2.95 ± 0.01. 

The two-dimensional ISG with Gaussian random in- 
teractions does not order above zero temperature^ Our 
data show that x(T) tends to 2.37 ± 0.06 at zero tem- 
perature, consistent with the value of 2.35 ± 0.02 ob- 
served in Ref. [J Both estimates suggest that E (O) 
is slightly more negative than the accurately measured 
6> DW = -0.28 ± O.Ol. 29 ! 30 ! 31 x{T) shows a shallow mini- 
mum atT« 0.4, unrelated to any ordering temperature. 

In the d = 2 GG, for which T c = 0, 24 x(T) increases 
steadily as T rises, from a value x(0.13) = 2.40 ± 0.02 
at the lowest measuring temperature T — 0.13 corre- 
sponding to #e(0) « —0.40 ± 0.02 in agreement with re- 
sults from Ref. 5J. Note also that the above estimate of 
x(T — > T c = 0) agrees with the zero-T c scaling of Eq. J3J). 

We can review the estimates for 6>e(0) in the GG. 
In dimension 2, the estimated E (O) = -0.40 ± 0.02 
could be consistent with either the "best twist" Obt, 
which is — 0.39±0.03j£*i 4 i or with the "random twist" 0rt, 
— 0.45 ±0. 0151^ For the three-dimensional GG, however, 
the low-temperature limit # E (0) = 0.010 ± 0.012 is in 
agreement within the error bars with the 9b,t estimate 
by Akino and Kosterlitz, 0.05±0.05 (Ref.[ll|) but is com- 
pletely different from the domain- wall 6>bt(0) estimated 
to be 0.27±0.01i£i±i For the GG in dimension 4, our data 
do not go below T = 0.7; we estimate E (O) = 0.54±0.05. 
By extrapolation from dimensions 1 (Ref. 32 : ), 2 and 3, 
we estimate for the four-dimensional GG 0bt ~ 0.9±0.1. 
We conclude that in the GG 0e(O), defined unambigu- 
ously through the effect of periodic boundary conditions 
on the energy, can be very different from the domain- 
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wall 6>bt(0); it may well be possible to identify it with 
#rt(0) as defined by Akino and Kosterlitz«ii Regardless, 
the size effect provides an operational physical realiza- 
tion of a bona fide 6 quite distinct from the "best twist" 
domain- wall stiffness exponent #bt(0). The strong dif- 
ference between the domain- wall exponent #bt(0) and 
the periodic boundary conditions exponent #e(0) seen in 
d = 3 and implied in d = 4 is presumably related to the 
extra liberty that vector spins (as opposed to Ising spins) 
have to reorganize under external constraints. Any dif- 
ferences there may be between 6>e(0) and #dw(0) in ISG's 
certainly appear to be less spectacular than in the GG. 
However, in both dimension 3 and dimension 2, inde- 
pendent measurements of #e(0) are consistent with each 
other while the agreement with #dw(0) is poor. 

We find that for the models studied, the observed crit- 
ical exponent x(T c ) is systematically higher than its rig- 
orous scaling value d—l/v. A possible explanation might 
be the influence of a "lattice artifact" corrections^ so 
that, for instance, e(L) - = aL _(d_1/l/) + bL~ d . In 
the presence of this correction term the data can be rep- 
resented quite accurately by a single effective exponent 
x with a value between d—l/v and d if b is positive. 

In conclusion, we have presented numerical results on 
the scaling correction to the internal energy per spin as 
a function of system size and temperature in a variety 



of spin-glass models. The T = finite-size correction 
has been linked to domain walls, and scaling predicts 
x = d — 1/v at T c . By definition, the domain-wall stiff- 
ness exponent #dw(T) drops to zero at T c (see Refs. ITU 
and l2lh and one would expect from the domain-wall 
picture that 9e(T) should drop in a similar way. Based 
on standard scaling arguments, the effective stiffness ex- 
ponent 9~e{T = T c ) should be equal to 1/v. In prac- 
tice, for < T < T c there is a steady enhancement of 
the effective stiffness exponent in all finite- T c spin-glass 
systems studied, which can be interpreted as a gradual 
change from domain-wall to critical behavior. Neverthe- 
less, the temperature-dependence of 9e is incompatible 
with a domain- wall mechanism. However, the critical 
dip is never as deep as the scaling theory predicts - the 
observed effective 9e(T c ) is systematically smaller than 
1/v. A rigorous relationship between the minimum of 
x(T), Tnim, and T c could provide a method to determine 
T c from purely energetic measurements. 
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